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We study the interaction between dark sectors by considering the momentum transfer caused by 
the dark matter scattering elastically within the dark energy fluid. Describing the dark scattering 
analogy to the Thomson scattering which couples baryons and photons, we examine the impact of 
the dark scattering in CMB observations. Performing global fltting with the latest observational 
data, we find that for a dark energy equation of state w < —1, the CMB gives tight constraints 
on dark matter-dark energy elastic scattering. Assuming a dark matter particle of proton mass, we 
derive an elastic scattering cross section of ao < 3.295 x W^^'^ar where ctt is the cross section of 
Thomson scattering. For u; > — 1, however, the constraints are poor. For w = —1, an can formally 
take any value. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 

There has been convincing evidence indicating that 
our universe is expanding at an increasing rate and it is 
mainly composed of dark matter (DM) and dark energy 
(DE) at the present moment. The simplest explanation 
of DE is the cosmological constant, that fits observations 
very well. However, the observed value of the cosmolog- 
ical constant falls far below the value predicted by any 
sensible quantum field theory. This is the so called cos- 
mological constant problem. Furthermore using the cos- 
mological constant to explain the DE unavoidably leads 
to the coincidence problem, namely, why are the vacuum 
and matter energy densities of precisely the same order 
today a? 

Considering that DE and DM contribute significant 
fractions of the contents of the universe, it is natural, in 
the framework of field theory, to consider the interaction 
between them. The covariant description of the interac- 
tion between DE and DM is given by Q Vf^Tl^^-^ = Q^^^, 
where Q'^x^ is a four vector governing the energy momen- 
tum transfer between DE and DM. The subindex A refers 
to DM or DE component respectively. For the whole 
system, the energy and momentum are conserved, and 
the transfer vector satisfies J2\ Q{x} = 0. The four vec- 
tor Q'^x) describing the interaction between DE and DM 
can be phenomenologically decomposed into two parts 
with respect to a given observer A' with four velocity 
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namely, Q'^^-^ 



F, 



(A) I 



where 



Q{X)\x' = —U(x')i^Q'(x) the energy transfer rate of the 
A component observed by the observer situated on the 
A' component. F(x)\^' ~ correspond- 
ing momentum transfer between the two components ob- 
served by the observer in the A' frame. 

The possibility that DE and DM interact with energy 
exchange between them has been studied extensively [J]- 



[3l| . It has been shown that the energy transfer between 
DE and DM can provide a mechanism to alleviate the co- 
incidence problem j3|-[lol|. Complementary observational 
signatures of the energy exchange between DE and DM 
have been obtained from the cosmic expansion history by 
using the WMAP, SNIa, BAG and SDSS data etcQiia]- 
[3] and the growth of cosmic structure p5j|-[3l|. 

Besides the energy transfer between DE and DM, their 
interaction may also impart a transfer of momentum. 
Considering the extremely low DE density, and the non- 
relativistic velocities of DM motions, elastic scattering 
appears naturally between the DM and DE fluid. We do 
not need to restrict to a particular DE model since the 
macroscopic physics is independent of the microphysics 
in the scattering process. The DM scattering elastically 
within the DE fluid results in the momentum transfer 
between them. The implications of the elastic scatter- 
ing between DE and DM was explored in [s^]- It was 
found that the growth of structure was suppressed by a 
drag term arising from elastic scattering between dark 
sectors. It is of great interest to extend the study in [sl] 
to investigate the cosmological signal of the elastic scat- 
tering between the DM and DE fluid. In this work we 
will explore the implications of the momentum transfer 
between dark sectors in the cosmic background radiation 
and constrain the scattering cross section between DE 
and DM by using the WMAP observation and the joint 
analysis of CMB, SNIa and BAG data. 



II. PERTURBATION FORMALISM WITH 
DARK SCATTERING 



In this section, we construct the first order perturba- 
tion theory when there is elastic scattering between the 
DM and the DE fiuid. The perturbed space-time at first 
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order reads 

ds'^ =a^[-(l + 2i:)dT^ + idiBdrdx' 

+ (1 + 2(t))6ijdx^dx^ + DijEdx'dx^] 



(1) 



where a is the cosmic scale factor, il), (p, B, E are the 
scalar metric perturbations and Dij = {didj — -1(5^ V^). 
Choosing the Newtonian gauge, we take B = E = 0. In 
the background the energy conservations of DM and DE 
are described by 



+ SHpc = aQc, 
p'd + 3?^(1 + w)pd = aQrf, 



(2) 
(3) 



where the subscripts 'c' and 'd' refer to DM and DE re- 
spectively. A prime denotes the derivative with respect 
to the conformal time r. w represents the dark energy 
equation of state and "H = ^ is the Hubble parameter. 
Here Qc and Qd are the energy exchange between DE and 
DM observed in background. If there is no energy trans- 
fer between the dark sectors, Qc = Qd = 0, the energy 
densities of DM and DE are separately conserved. 

In the Fourier space, the covariant form of the per- 
turbed energy-momentum transfer between DE and DM 
can be expressed as [IF 
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5'^ + 3H(^ - wx)5x + + wx)kvx 
opx 

= -3(1 + wx)<l>' + (2V. - Sx)^ + 



v'x + T-L{1 - 3wa)wa - -r— — t-^^x 



l + wx Spx 



1 + Wx 



VX + kil; ■ 



Px [i + wx)px 



a oQpx 



(4) 
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where Sx = dpx/px and vx is the peculiar velocity. 6Qpx 
is the potential of the spatial part of the perturbed cou- 
pling vector, 6Q\. As discussed above, 6Qpx can be de- 
composed into two parts, SQpx = Qx\x'Vx' + ^a|a'- vx' 
is the velocity of the observer situated in the A' com- 
ponent and Fx\x' is the external non-gravitational force 
density on A component observed by the observer sitting 
in the A' component. In the following discussion, we will 
assume that the energy transfer vanishes in the back- 
ground and only concentrate on the momentum trans- 
fer between dark sectors. Therefore, in the background 
frame, Qx\{back) = 0, and only Fx\{back) contributes to 
the energy momentum transfer. Fx\(^back) is simplified as 
Fa below. 

In the absence of a fundamental theory, we can only 
conceive a phenomenological description of DM scatter- 
ing elastically within the DE fluid. Before we do so, we 
recall the description of the coupling between baryon and 
photon fluid. The perturbation equations of baryon read 



S'^ + iClndb + 30' = -kvb 



(6) 



v'u = — "Hwfc + fcCf(5b + fc^/j + anf.(TT-r-^{v-y — t'b), (7) 
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where the subscript represents baryon and '7' stands 
for photon. is the sound speed. The last term in ([7]) 
arises from Thomson scattering. Comparing ([6]) and ([7]) 
with dU) and ©, we found that Q° = 0, 5Ql = and 



Fb = 7ieCTT(l + w^)p^{v^ - Vb)/a. 



(8) 



Here ctt is the cross section between baryon and photon. 
Using the analogy, we can choose the interaction term of 
DM due to the elastic scattering in the DE fluid 



Fc = -{1 + w)pdnc(JD{vd - Vc), 
a 



(9) 



where ajy represents the unknown cross section of the 
elastic scattering between DM and DE, Uc = pc/mc is 
the number density of DM particles. Conservation of 
momentum leads to a similar term arising for DE. It is in- 
teresting that {1+w) appears in the interaction term. We 
can see that this is a general requirement for momentum 
transfer. Given A refers to dark energy va Fd — —Fc 
is divided by (I + w), which must vanishes when w ap- 
proaches — 1 in order to avoid singularities. However, this 
does not imply that momentum transfer has no impact 
on the dark sector perturbations for w ^ —I. Since we 
do not have a good estimation of the mass of DM par- 
ticles, it is convenient to define S = ctd/wc and express 
the momentum transfer as 



Fc = -{1 + w)pdPc^{vd - Vc). 
a 



(10) 



This is the ansatz of the dark scattering adopted in |32] ■ 
In the study of the scattering between baryon and dark 
fluid, the same interaction term for the dark fluid was 
adopted in ^33]. It would be fair to say that we do not 
know the microphysics in quantizing the DE so that we 
do not have the exact definition of the cross section of 



the elastic scattering between DM and DE. In [32| it was 
argued that the analysis of the macroscopic behavior are 
largely independent of the microphysics involved and the 
bound on the DM and DE cross section was first derived 
from the impact incurred on the growth of large scale 
structure. In this work we are going to investigate the 
signature of this phenomenological term on the DM and 
DE elastic scattering in the CMB observations. For the 
sake of simplicity, we will set S to be constant in the 
following. 

Inserting ^ into the general equations ^ and JS]), we 
get, in the Newtonian gauge, the evolution equations of 
the perturbations to DM and DE 

S'c = - kvc - 30', (11) 
v'c=- Uvc + kijj + a{l + w)pdT;{vd - Vc), (12) 



Si, 



Vd 



Spd 
Spd 



w)Sd - (1 + w)kvd 



-3(l + w)0', 
k 6pd 



(13) 



1 + w Spd 
+ apcT,{vc - Vd) 



w 

Sd - H(l - 3w)vd - — — Vd + kijj 
1 + w 



(14) 
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We choose the Newtonian gauge where the Thomson 
scattering is weh estabhshed. Now we rewrite the equa- 
tions in the gauge invariant quantities so that our calcu- 
lations below do not depend on the specific gauge choice. 
Constructing the gauge invariant quantities [lOj] 
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we obtain the gauge-invariant perturbation equations 



D' = 



-kVc , 

- HVc + fc* + a(l + w)pd^{Vd - Vc) 
-k{\ + w)Vd-m{Cl-w)Dd 

+ {3w' + 9H{C^ ~ w){l + w)}<i> , 
w' 

- ^—Vd - n{\ - ?,w)Vd + kcl- 

l + w 1 
+ m{Cl -Cl)Vd- 3A:C2$ + k^ 
+ ap,^[Vc - Vd) . 



Dd 



We have employed 

Cl5d - iCl 



Spd 
Pd 



B 



Pd 



(15) 
(16) 



(17) 



(18) 



(19) 
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where Cg is the effective sound speed of DE and 

is the adiabatic sound speed of DE. We assume Cj 
in our calculation below. The elastic scattering between 
DM and DE does not enter the equations governing the 
perturbations of energy densities, but affects the pertur- 
bations in the peculiar velocities. 



A. Stability of the perturbation 

It is known that DE perturbation is gravitationally 
unstable when its equation of state crosses —1 [34]- 36]. 
Such an instability happens regardless of whether there 
is interaction between DE and DM or not. To avoid 
this instability we restrict 1 -I- w to be either positive or 
negative in the following discussion. 

When there is energy transfer between DE and DM, 
it was observed that the stability of the perturbation de- 
pends on the form of the interaction and the equation of 
state of DE (1 111- [11. 



It is of interest to examine the stability of the pertur- 
bation when there is momentum transfer between DM 
and DE. We find that for the DE equation of state satis- 
fying l-l-ui > 0, there is no instability in the perturbation 
caused by the scattering between dark sectors. The dark 
scattering acts as a new drag term in the perturbation 
equations. In the early period, when apc dominated and 
apd was small, we observe that Vc evolved freely while Vd 
was forced to follow along with Vc- With the decrease of 
apc Vd gradually departed from Vc and decayed. In the 
late time, when a(l -|- w)pd became important, the dark 
scattering will drag Vc to Vd- The effect of the scattering 
between DE and DM on the behavior of the evolution of 
Vc, Vd can be clearly seen in FiglTal 

In the case when {1 + w) < 0, the sign of the last term 
in (fT6|) is flipped and the drag force turns into propulsion. 
This will cause the blow up of the peculiar velocities of 
DM and DE, which are shown in Figllbl To understand 
the blow up of Vc, we can simply look at 1/' ~ a{l + 
w)pd'S,{Vd — Vc)- Considering \Vd\ < \Vc\ in the late time, 
w < —1 will cause the exponential growth of Vc- In turn 
the exponential growth of Vc will also lead the blow up 
of Vd due to the last term in p8|) . The instability due to 
dark scattering is driven by the term a(l + w)pd^{Vd — 
Vc), which is suppressed when w is close to —1. This can 
be seen clearly in FigfTcl 



B. Tight coupling approximation 

Since we do not consider the energy transfer between 
DE and DM, the DM energy density pc scales as 
and the combined term apc^ scales as a~^- The over- 
whelmingly large value of apc^ in the early time becomes 
the bottleneck in the computation which severely limits 
the step size and the speed in the numerical calculation. 
It precludes a straightforward numerical integration at 
early times: tiny errors in the propagation of Vc and 
Vd will lead to serious errors to the force. This prob- 
lem is similar to that encountered when the baryon and 
photon fluids were str ong ly coupled via Thompson scat- 
tering at early times [33] • One resorted to the so called 
tight coupling approximation to improve the numerical 
computation in CMB. Here we will borrow this method 
and eliminate the term of order apc and apd from the 
evolution equation. 

Subtracting (fT8| from (fT6|) and doing the rearrange- 
ment, we have 



a[(l 



w)pd 



1 



- Pc]SS' = -S' 
+ Ut\ Vd + kCl 



HVc 
Dd 



1 



3fcCf$, 



(20) 



where 5 = - K and T = -(1 - 3w) + ?,{Cl - CD- 
Taking derivative with respect to the conformal time, we 
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FIG. 1: The evohitions of pecuhar velocities for 
fc — O.lMpc^"'^. S is in the unit of aT/iTip where ctt is 
the cross section of Thomson scattering and nip is 
proton mass. In (a) and (b), thick hnes stand for Vc 
and thin hnes represent Vd- In (a) we choose w — —0.8 
and in (b) w = —1.5. In (c) we show the behavior of 
pecuhar velocity evolution of Vc when w approaches — 1 
from below when we fix S = 1 x IO^^ctt/™?)- 



«'[(! + w)pd + Pc]S5 + a[{(l + w)pa}' + p^^S 
+ a[{l+w)pd + Pc]^S' 

= -S" + H'Vc + HV/. 
+ H'T + UT'] Vd 
-kC? 



l + w 



1 + w (1 + wf 

nAv^ 



(21) 
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The derivatives appearing above can be written as 
p'c — —SHpc , p'd = — 3H(1 + w)pd , a' = oH . 

Inserting these expressions together with ([T5|) and ([1 
in (|2T]). it becomes 



a[(l + w)pd + Pc]SS" 

^^S" -[[l + w)pd-pc]-^<- 

l + w 

+{3(1 + wfpd - (T - 2)p Ja'HSS' + WV^ 



-Wv, 



l + w 



I + IU (1 + w)2 
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-3fcCf 
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Dd 



(22) 



+'HT]kCl- 

l + w J ''l + w 



$ - ikClUT^ - 3fcCf $' 
+ H{l + T)\k-^. 



Inserting (1^^ into ([^0]) and in turn using in ([T5)) and 
(IT51) . we obtain the evolution equations for Vc and Vd, 
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l + w 
I + w + r 



S' - 3/cC> 



l + w 
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+ -HrVd + fcC',^ 



■HT I Vd + kCi 
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(23) 



l + w 
- 3fcC2$ 



1 + w + r 
w' 



l + w 



UTWd + kC, 



l + w 

S' - ikCl<^ + nVc (24) 
Dd 



l + w 



where r = Pc/ Pd- There is no longer any term propor- 
tional to apc or apd in the above formulae. Up to now, we 
have not yet made any approximation. To finish the tight 
coupling approximation, we drop the term S" in p2|) as 
has been done in [37[. The second derivative of S is the 
only unknown term during the numerical integration of 
perturbations and can be neglected in the expression of 
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S' when DE and DM are strongly coupled. In the tight 
coupling approximation we adopt here, considering the 
DE density to be sufRciently low at early times, Vc can 
evolve along (ITS)) while for Vd we need to employ the ap- 
proximated equation (j24p . in which S' is evaluated using 
([22|l and S" therein is neglected. This strategy can lead 
to the efhciency in computation. 



III. IMPACT OF THE DARK SCATTERING ON 
THE CMB POWER SPECTRUM 

With the perturbation equations at hand, we can study 
the influence of dark scattering on the CMB radiation. 
We have modified the public available code CMBEASY 
to account for the elastic scattering between DE and 
DM. In our theoretical computation, we allow the equa- 
tion of state of DE to be constant and time-dependent 
in the range w > —1 and w < —1 respectively. To see 
the effect of the scattering in CMB, a wide range of the 
strength of dark scattering is explored for the selected 
DE equation of state. 

When w is very close to —1, the momentum trans- 
fer due to the elastic scattering between DM and DE 
does not lead to a clearly discernible change on the CMB 
power spectrum. The dark scattering effect can obviously 
show up once the equation of state of DE deviates more 
from —1. For a constant equation of state of DE, the dark 
scattering impact on CMB can be seen in Figl2l When 
the constant equation of state of DE w > —1, we find 
that with the increase of the strength of the scattering, 
the small I CMB spectrum is enhanced. This property 
does not hold for constant DE equation of state satisfy- 
ing w < —I. In Figl2b| it is shown that the low I CMB 
spectrum will first be suppressed and then enhanced with 
the increase of the strength of the dark scattering. When 
the scattering is strong enough, the instability we encoun- 
tered in the perturbation will cause the blow up in the 
low I CMB spectrum. In Fig|3]we show the behavior of 
the gravitational potential when there is dark scattering 
between DE and DM. The gravitational potential is the 
source term for the low I CMB spectrum. For a constant 
equation of state of DE w > —1, a larger value of the 
scattering results in a further change of the gravitational 
potential. When the constant DE equation of state sat- 
isfies w < —1, there is a sharp increase in the potential 
as shown in Fig jSb] when the scattering is strong enough, 
which accounts for the blow up in the large scale CMB 
power. Since observations suggest the deficit of large 
scale CMB power, thus we can rule out the too strong 
scattering between DE and DM. Except for the imprints 
in the low I CMB spectrum, different from the energy 
transfer between dark sectors we do not see any 

change on the acoustic peaks due to the dark scattering 
when we choose the DE equation of state to be constant. 

Besides the constant equation of state of DE, we also 
examine the effect of the dark scattering on CMB when 
w is time-dependent. We choose the Chevallier-Polarski- 
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FIG. 2: CMB angular power spectrum when there is 
dark scattering between DE and DM. S is in the unit of 
ctt/tip. In (a) we choose constant DE equation of state 

w = —0.8. In (b) the constant DE equation of state is 
fixed as w = —1.5. 



Linder (CPL) parametrization|39| to describe the time 
varying equation of state of DE, where w is expressed 
as w{a) — au'o + (1 — a)we- In the early time, a ^ 1, 
w ~ We; at present, a ~ 1, w ~ wq. We restrict the 
equation of state of DE to be in the range either w > — 1 
or w < —1 to avoid the singularity at w = — 1 in the 
perturbations. 

When the DE equation of state is time-dependent and 
obeys w > —1, on the low I CMB spectrum, we observe 
a similar effect of the dark scattering to the case for con- 
stant DE equation of state satisfying w > —1, see FigHal 
In addition, the signature of the momentum transfer due 
to the elastic scattering between dark sectors appears in 
the first acoustic peak as well. The change of the acous- 
tic peak mainly appears for a large value of We, however 
this effect is washed out when is small. 
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FIG. 3: The evolution of gravitational potentials when 
k — O.lMpc"^. In (a) we choose w = —0.8 and in (b) 
we select w = —1.5. 
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FIG. 4: CMB angular power spectrum of models for 
dynamic DE equation of state in the form 
w{a) — woa + We{l — a). In (a) we choose 
Wo — —0.999, We = —0.2 and in (b) we select 
Wo = -1.001, We = -1.5. 



When the CPL form of the DE equation of state is 
in the range w < —1, the instability can also appear 
as the case of constant DE equation of state. In the 
CPL parametrization, We relates to the DE equation of 
state at the early time, while wo relates more to the late 
time DE equation of state. From Figl5]we observed that 
either wo or We deviating significantly from —1 will lead 
to the blow up in the CMB power spectrum at small 
I. The blow up is caused by the quick growth of the 
velocity when the DE equation of state is much smaller 
than —1. Choosing appropriate DE equation of state, 
we can obtain similar dependence of the CMB spectrum 
on the scattering between dark sectors observed for the 
constant DE equation of state as shown in FiglH The 
strength of the dark scattering is constrained to suppress 
the CMB spectrum at small I. 



IV. FITTING RESULTS 

In this section wc confront our models with observa- 
tional data by implementing joint likelihood analysis. We 
take the parameter space as 

P = {h, UJb, i^m, T, \n[10^" As], Us, Wo, We, ^) 

where h is the bubble constant, ujh — fli,h^ and u!„i — 
flmh^ are baryon and matter abundances, r is the opti- 
cal depth to last scattering surface. As is the amplitude 
of the primordial curvature perturbation, Us is the scalar 
spectral index. E here marks the strength of the dark 
scattering, which is in the form E = au/nic, where ao is 
the cross section of Thomson scattering between DM and 
DE fluid, rric should be the mass of DM particles. Since 
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FIG. 5: CMB angular power spectrum of models with 
different time-dependant DE equation of state. The DE 

EoS is in the form w{a) = woa + ?iie(l — a). 
E = 1 X 10^^(TT/mp. In (a) we fix We = —1.5 and in (b) 
we select wo = — f .001. 
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FIG. 6: The ID likelihood of the coupling strength of 
dark scattering, S. Thin solid lines correspond to fitting 
results using WMAP data only. Thick lines correspond 
to the result of fitting combined data. In (a) we choose 
w > —1 and in (b) we deal with w < —1. 



there is no good approach for computing the DM parti- 
cle mass, we fit S from observational data, instead of the 
cross section aj). We will choose the CPL parametriza- 
tion for the DE equation of state and consider the spa- 
tially flat universe with = 0. To avoid the singularity 
at w = —1 in the perturbation, we concentrate on the 
range of the DE equation of state satisfying w > —1 
and w < —1, respectively. Our numerical computation is 
based on the CMBEASY code [si]. We will first fit the 
CMB anisotropy data from the seven-year Wilkinson Mi- 
crowave Anisotropy Probe(WMAP). The fitting results 
are listed in TableU) The coupling strength E is in the 
unit of (Tt /mp where ctt is the cross section of Thomson 
scattering and nip is proton mass. The Id marginalized 
likelihood for the strength of the dark scattering is shown 
in the thin solid lines in FiglHlwhen we only use the CMB 



anisotropy data from the seven-year WMAP observation. 

In order to get tighter constraint on the parameters, we 
employ BAO and SNIa data sets in addition to WMAP. 
The BAO distance measurements [i^l which are obtained 
from analyzing clusters of galaxies and test a different re- 
gion in the sky as compared to CMB. BAO measurements 
provide a robust constraint on the distance ratio 

d,^r4zd)/D4z) (25) 

where Dy{z) = [{1 + z)'^D\z/ H{z)]^/^ is the effective dis- 
tance [111, Da is the angular diameter distance, and H{z) 
is the Hubble parameter. rs{zd) is the comoving sound 
horizon at the baryon drag epoch where the baryons de- 
coupled from photons. We numerically find Zd using the 
condition t/R = 1,R = as defined in i42j. The 
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TABLE I: The cosmological parameters from the global fitting. 





WMAP 


WMAP- 


hBAO+SN 


Parameter 


w > —1 w < —1 


w > -1 


TO < -1 



h 

T 

n 

ln[10"'A, 
wo 

We. 
E 



n riT?! +0 0008 

U.UZ01_Q 0007 

0.7241°°^,^ 
0.089lg°}^ 

00+0 024 
Q nc7+0.036 

< -0.948 

< -0.943 

N/A 



n 1Q4+0.005 
U.iOy:_o.o05 

U.UZZU_o.oo06 

u. roz_o 031 
083+° °" 

U.UOO_o 013 

n qc^c^+0.014 
U.»00_o 014 

q n7'?+0 034 
■J-"' '5-0.034 

> -1.260 

> -1.370 

< 3.295 X 10~"' 



0-13llo: 
0.0228to 
0.706^°: 



0.0841 
0.982^ 



3.070t^ 



004 
004 

,0007 
,0006 

013 
013 

).015 
).014 

).019 
).017 

,035 
,033 



< -0.950 

< -0.960 
N/A 



136+0 004 
U.±OU_o.o04 

0.0220t°:S 

0.698tHii 
082+° °" 

U.UOZ_o 013 

n Qc:tS+0 013 
U.»00_o 013 

Q n7fi+° °33 

> -1.031 

> -1.077 

< 1.665 X 10"® 



X%Ao calculated as [4C 



-tobs 



(26) 



where d = (4=0.2,4=0.35)^, d = (0.1905,0.1097)^ 
and the inverse of covariance matrix reads 14011 



30124 -17227 
-17227 86977 



Furthermore, we add the BAO A parameter (4]J 



(27) 



A 



i;(0.35)i/3 [0.35 
0.469(ns/0.98) 



0.35 



dz 







-0.35 



E{z) 
±0.017 



2/3 



(28) 



where E(z) = and Us are the scalar spectral index. 

We also use the compilation of 397 Constitution sam- 
ples from supernovae survey [i^. We compute 



XsN 



E 



and marginalize the nuisance parameter. 
We implement the joint likelihood analysis. 



X 



2 ,2,2 

XWAIAP + XsN + XbAO- 



(29) 



(30) 



From FiglSal we see that CMB anisotropy probes can 
not put tight constraints on the intensity of dark scatter- 
ing when w > —I. One reason to understand this is that 
the dark scattering mainly influences the small I CMB 
power spectrum through the late ISW effect, while the 
CMB data at small I are poor and lacking. Further, from 
the theoretical study we learnt that the elastic scattering 
between DM and DE is suppressed when the equation of 
state of DE is close to —1. From the fitting results for 
the equation of state of DE, it is really close to —1 in the 
recent epoch, w ~ — 1 is further constrained by adding 



BAO or SNIa data. When the DE equation of state is 
in the range w < —1, we see from Fig|6bl that the dark 
scattering is tightly constrained. This is because when 
w < — 1 the small I CMB spectrum is enhanced when 
the scattering is strengthened and blows up if the scat- 
tering is strong enough. The observed suppression of the 
CMB data at low / strongly puts the upper bound on the 
strength of the scattering. For this reason, a tight con- 
straint is put on the dark scattering when the DE equa- 
tion of state obeys w < —1. Including the BAO or SNIa 
data, the equation of state of DE is constrained further 
close to —1, the cross section (Jd can then be arbitrarily 
large. All we can really constrain is E' = (1 -|- w)(Jd /nric. 
There is a degeneracy between {1+w),(Jd and rric- When 
w ^ — 1 is constrained at present, the constraints on ajj 
will be poor. 

If we assume rtic — nip, where vrip is the mass of proton, 
from the fitted value of S in Table. 1, we can estimate 
UD < 3.295 X 10-i°CTT (for WMAP alone) and (Jd < 
1.665 X 10~^ctt (from the combined constraint). 

To conclude, when w > ~1, the effect of elastic scat- 
tering between dark sectors has little influence on the 
CMB power spectrum, which is different from the en- 
ergy transfer discussed in our previous papers. Thus the 
CMB observation cannot put tight constraint on the dark 
scattering when DE equation of state w > —1. When 
w < —1, in the constrained ranges of wq and We, big 
strength of the dark scattering will lead to the sharp blow 
up of CMB spectrum at small I. This effect is similar to 
the energy transfer discussed in Q [I^ H ■ Using the pre- 
cise CMB data, the strength of the dark scattering can 
be constrained tightly. 



V. CONCLUSION AND DISCUSSION 

In this paper we generalized the study of the inter- 
action between DE and DM. Instead of considering the 
energy exchange between DE and DM in the literatures. 
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in this work we investigated the momentum transfer by 
exploring the elastic scattering between dark sectors. 
Adopting the phenomenological description of the DM 
elastically scattering within the DE fluid '32'], which is 
analogous to the Thomson scattering for baryons and 
photons, we examined the impact of the dark scattering 
on CMB observations. We found that in addition to the 
growth of structure [s^l, the elastic scattering between 
the DM and DE fluid can also leave cosmological sig- 
nals in CMB. The imprint of the dark scattering in CMB 
highly depends on the equation of state of DE. When 
the DE equation of state w < —1, we can constrain the 
dark scattering well from CMB observations. However 
when w > —1, the cosmological signal of the scattering 
is weak. The impact of the dark scattering is suppressed 
when w —1. 

From our result, we learnt that the CMB probes are 
not sufficient to constrain the dark scattering. It would 



be interesting to carry out the complementary constraints 
by combining the CMB test with the large scale struc- 
ture tests as employed in the study of the energy transfer 
between dark sectors [25'|-'3l]. Progress in this direction 
will be reported in the future. 
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